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1.  INTRODUCTION 


The  technical  objective  of  this  research  project  is  to  develop  flight  dynamics  analysis  methodol¬ 
ogy  for  combat  damaged  fixed-wing  aircraft  to  support  ballistic  vulnerability  assessments.  With 
combat  damage,  wing  structural  capability  (stiffness/strength)  can  be  affected.  When  critical 
structural  members  are  weakened  or  destroyed,  aerodynamic  loadings  may  cause  the  structure  to 
fail.  In  this  study,  removal  of  wing  section  is  considered. 

Removal  of  wing  structure  caused  weight  unbalance  and  the  loss  of  the  aerodynamic  lifting 
surface.  Due  to  the  loss  of  a  wing  mass,  a  change  of  c.g.  (center  of  gravity)  occurs.  The 
aerodynamic  configurations  of  damaged  aircraft  generally  involve  unusual  shapes,  asymmetrical 
configurations,  and  possibly  large  angles  of  attack  and/or  sideslip  for  trimmed  flight  Previously, 
standard  configuration  codes  have  been  used  to  calculate  aerodynamic  stability  derivatives.  For 
example,  the  USAF/DATCOM  and  AAA  (Advanced  Aircraft  Analysis)  codes  have  been  used  for  this 
purpose.  Input  for  these  codes  comes  from  the  geometry  and  the  physical  properties  of  the  aircraft 
configuration  under  study.  The  Air  Systems  Branch  of  US  Army  Research  Laboratory  has  also 
investigated  the  capability  to  calculate  aerodynamic  control  and  stability  derivatives  and  trim 
values  for  damaged  configurations  using  the  MSC/NASTRAN  aeroelastic  code,  a  system  of 
computer  programs  for  modeling  aircraft  aerodynamics  and  structures,  and  for  analyzing  aircraft 
stability  and  loads. 

The  finite  element  structural  model  and  aerodynamic  models  of  the  aircraft  are  developed 
individually  and  splined  for  the  interpolation  of  deflections.  Splines  for  both  the  lines  and  surfaces 
are  used  to  generate  the  transformation  matrix  from  the  structural  grid  point  deflections  to  the 
aerodynamic  grid  point  deflections,  where  local  streamlines  are  also  computed.  In  this  study,  all 
wing  skins,  stabilator  skins,  vertical  wing  skins,  and  fuselage  skins  are  considered  as  quadrilateral 
plate  elements,  and  the  surface  spline  interpolation  method  is  used.  Of  particular  interest  is  the 
wing-body  interference  option  used  in  this  report  as  compared  to  results  in  [1],  where  this  option  is 
not  used. 

In  this  study,  four  different  subsonic  Mach  numbers  are  considered  to  analyze  the  static 
stabilities,  aerodynamic  influence  coefficients,  and  trim  values  in  level  flight.  Static  longitudinal  and 
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lateral-directional  stabilities  are  calculated  and  discussed  for  the  undamaged  and  damaged  wing 
configurations.  The  aircraft  body  is  treated  as  elastic  and  unrestrained  (free-free  boundary  condi¬ 
tion).  Aerodynamic  influence  coefficients  are  used  to  calculate  the  aerodynamic  pressure, 
aerodynamic  forces,  and  moments  at  subsonic  speeds. 

2.  STATIC  AEROELASTICITY 

The  static  aeroelastic  analysis  involves  the  responses  of  a  flexible  structure  to  aerodynamic 
loading  and  yields  the  static  response,  stability  and  control  derivatives,  trim  variables,  air  loads,  and 
stresses  and  strains  of  the  wing  structure.  All  these  variables  are  obtained  for  an  unrestrained 
(free-flying)  aircraft  structure.  The  MSC/NASTRAN  aeroelastic  solution  uses  an  aerodynamic 
influence  coefficient  matrix  that  is  generated  from  data  describing  the  geometry  of  the  aerodynamic 
finite  elements.  The  aerodynamic  influence  coefficients  from  the  Doublet-Lattice  method  are  used 
for  the  calculation  of  aerodynamic  quantities  in  subsonic  flow.  This  is  a  linearized  aerodynamic 
potential  flow  theory  which  is  presented  in  [  2-4].  The  Doublet-Lattice  method  with  body  interference 
is  only  used  for  the  subsonic  ranges.  The  details  of  the  static  aeroelastic  analysis  is  explained  in 
Appendices  A-C. 

3.  GENERAL  APPROACH  TO  THE  PROBLEMS 

The  MSC/NASTRAN  aeroelastic  solution  requires  a  well-designed  aircraft  structural  model 
with  stiffness  (for  static  analysis)  and  mass  balance  (for  dynamic  analysis),  and  an  aerodynamic 
model  to  obtain  the  best  result.  The  reference  geometry  for  the  aerodynamic  and  structural  input 
data  is  also  required. 

To  analyze  the  aerodynamic  effects  of  a  damaged-wing  aircraft,  five  different  finite  element 
wing  models  are  examined.  The  first  is  the  baseline  (undamaged)  and  the  second  through  fifth  are 
for  damaged  wings.  For  the  damaged  cases,  the  right  wing  is  examined.  A  full-span  aerodynamic 
and  structural  model  of  the  wing  accounts  for  the  asymmetrical  configuration  of  the  aircraft. 
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3.1  Description  of  Aircraft  Model 


The  T-38  aircraft  construction  details  used  to  model  the  various  structural  assemblies  are 
described  in  [5]  and  [6].  See  Figures  1  and  2a-c.  The  cantilevered  wing  features  an  aspect  ratio  of 
3.75,  no  dihedral  or  angle  of  incidence  relative  to  the  fuselage,  and  a  sweepback  angle  of  25°  at 
quarter  chord.  The  wing  section  has  a  NACA  65A004-8  (modified)  airfoil.  The  mean  aerodynamic 
chord  is  c  =  93  in.,  the  chosen  reference  chord  is  C  -  80.26  in.,  and  the  thickness/chord  ratio  is  4.8 
%.  The  reference  wing  span  is  303.0  in.,  and  the  reference  area  of  the  full-span  wing  is  S  =  2431 9.1 
in2.  The  horizontal  stabilator  has  an  aspect  ratio  of  2.85  with  a  mean  aerodynamic  chord  of  43.63 
in.  The  reference  area  of  the  horizontal  stabilator  is  S  =  4610  in2.  The  fuselage  is  composed  of  an 
aluminum  semi  monocoque  basic  structure  with  steel,  magnesium,  and  titanium  sub-members.  The 
fuselage  length  is  L<  *  582.1  in.  The  weight  of  the  wing  is  1,170  lb.,  and  the  total  weight  of  the 
aircraft  is  11,500  lb.  The  moment  of  inertia  of  aircraft  is  lx*  =  6.86  x  107(lb-in2),  lyy  =  1.297  x  10® 
(lb-in2),  and  =  1.343  x  10®  (lb-in2). 


3.2  Structural  Wing  Model 

Five  finite  element  wing  models  are  created  in  this  report.  Four  damaged  wing  models  are 
modified  from  the  model  baseline  (undamaged  aircraft).  These  are  created  by  removing  various 
regions  (portions)  of  the  aerodynamic  surfaces.  Damage  case  1  is  created  by  removing  the  internal 
structure  and  top  and  bottom  skins  of  wing  tip  portions  between  w.s.  125.0  andw.s.  151.5.  Damage 
case  2  represents  more  internal  damage;  additional  spar  and  rib  webs/chords  and  skins  are 
removed  between  w.s.  125.0  and  w.s.  111.0.  Damage  case  3  has  additional  spar  and  rib 
webs/chords  and  skins  removed  between  w.s.  111.0  and  w.s.  76.7.  Damage  case  4  has  the  aileron 
control  surface  alone  removed.  These  four  damage  cases  are  shown  in  Figures  2d-g.  The  structural 
elements  removed  from  the  baseline  finite  element  model  to  represent  damage  are  listed  in 
Table  1.  Wing  damage  reduces  both  the  wing  area  and  the  lateral  distance  from  the  aircraft 
centerline  to  the  wing  centerline  of  pressure. 

To  analyze  the  maneuvering  flight  of  an  aircraft  with  a  damaged  wing,  the  wing  structural 
models  (damaged  and  undamaged)  are  added  to  the  fuselage  structural  model,  and  to  the  tail  wing 
structural  model  group. 


Table  1.  Structural  Element  and  Removed  Numbers  for  Each  Model 


Element  Numbers 

'"'\MODEL 

Baseline 

Damaae  1 

Damaae  2 

Damaae  3 

Damaae  4 

ROD 

343 

343 

343 

343 

343 

BEAM 

276 

276 

256 

214 

276 

SHEAR 

136 

136 

125 

105 

136 

CONROD 

89 

89 

81 

67 

89 

QUAD4 

522 

516 

510 

492 

520 

BAR 

30 

30 

30 

30 

30 

TRIA3 

104 

104 

104 

104 

104 

mm 


3.3  Wing  Aerodynamic  Model 


The  wing  aerodynamic  model  is  developed  using  MSC/XL,  initially  for  the  aerodynamic 
analysis  and  finally  (interpolated  by  the  surface  SPLINE  option  card  in  MSC/NASTRAN)  for  the 
aeroelastic  analysis.  In  these  analyses,  five  aerodynamic  box  patterns  are  used;  the  description  of 
input  data  and  the  aerodynamic  modeling  technique  is  explained  in  Appendix  A. 

Aerodynamic  models  for  undamaged  and  damaged  wing  cases  were  created  in  the  same 
manner  as  the  structural  wing  models.  As  shown  in  Figure  3,  five  cases  were  examined.  Region  1 
is  defined  within  w.s.  0.0  (aircraft  centerline)  and  w.s.  76.7,  and  within  spar  stations  0.0%  and 
100.0%.  Similarly,  region  2  is  within  w.s.  76.7  and  w.s.  111.0,  and  spar  stations  0.0%  and  66.6  %. 
Region  3  is  within  w.s.  76.7  and  w.s.  111.0,  and  spar  stations  66.6%  and  1 00.0%.  Region  4  is  within 
w.s.  111.0  and  w.s.  125.0,  and  spar  stations  0.0%  and  100.0  %.  Region  5  is  within  w.s.  125.0  and 
w.s.  151 .5,  and  spar  stations  0.0  %  and  100.0  %. 

Region  5  is  removed  for  damage  case  1.  Regions  4  and  5  are  removed  for  damage  case  2. 
The  wing  tip  station  is  changed  from  151 .5  to  125.0  for  case  1  and  to  111.0  for  case  2.  Regions  2, 
3,  4,  and  5  are  removed  for  damage  case  3.  Region  3  is  removed  for  damage  case  4.  For  the 
damage  case  4,  the  wing  tip  station  is  changed  from  1 51 .5  to  76.7.  Aerodynamic  boxes  are  located 
on  the  remaining  wing  and  tail  surfaces.  The  corresponding  aerodynamic  models  are  shown  in 
Figures  4a-g. 

3.4  Wing-Body  Interference 

Since  this  study  is  concerned  with  lifting  surfaces,  the  fuselage  is  represented  as  a  beam  and 
the  aerodynamic  effects  of  the  fuselage  are  neglected  (i.e.,  it  is  assumed  that  there  are  no 
aerodynamic  forces  acting  on  the  fuselage). 

The  fuselage  body  is  further  idealized  as  slender  with  interference  elements.  The  primary 
purpose  of  the  slender  body  elements  is  to  account  for  the  forces  arising  from  the  motion  of  the 
fuselage,  whereas  the  interference  elements  are  used  to  account  for  interference  among  all  bodies 
and  panels  in  the  same  group.  This  is  done  by  providing  a  surface  through  which  the  boundary 
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Figure  3.  Designation  of  removed  region  of  wing  structure. 


condition  of  no  flow  is  imposed.  Bodies  are  further  classified  as  to  the  type  of  motion  allowed.  In 
the  aerodynamic  coordinate  system,  y  and  z  are  perpendicular  to  the  flow.  In  general,  bodies  may 
move  in  both  the  y-(lateral)  and  z-(vertical)  directions.  Frequently,  a  body  (i.e.,  a  fuselage)  lies  on 
a  plane  of  symmetry  and  only  z-(or  y-)  motion  is  allowed.  Thus,  this  model  contains  zy-bodies.  One 
or  two  planes  of  symmetry  or  antisymmetry  may  be  specified.  Figure  5  shows  an  idealization  for 
determining  the  effect  of  the  wing  on  itself  and  other  parts  of  the  aircraft.  More  details  of  wing-body 
(fuselage)  interference  used  in  this  study  are  described  in  [7]. 

3.5  Control  Surfaces 

The  aircraft  control  surface  deflections  are  defined  by  the  conventional  symbols  8,,  8a,  and  Sr 
for  pitch,  roll,  and  yaw  control,  respectively.  The  pitch  control  deflection  5b  is  defined  as  a  symmetric 
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Figure  4a. 


Figure  4c. 


Figure  5.  Idealized  interference  wing  and  body. 


equal  deflection  of  the  left  and  right  horizontal  tail  surface  (stabilator),  while  the  roll  control  deflection 
5a  is  defined  as  an  antisymmetric  deflection  of  the  ailerons.  Yaw  control  deflection  5r  is  defined  by 
rotation  of  rudder. 

Positive  aileron  deflection  is  defined  as  left  aileron  trailing  edge  down,  right  aileron  trailing  edge 
up,  and  produces  a  positive  rolling  moment  Since  the  left  and  right  stabilators  rotate  together(one- 
piece  all  moving  tail)  the  positive  rotation  is  defined  as  the  leading  edge  down  for  both  sides  and 
produces  a  nose-up  control  moment  A  trailing  edge  left  rudder  rotation  gives  a  positive  side  force 
from  the  rudder.  The  control  deflection  5a(aileron),  ^(stabilator),  and  5r{ rudder)  are  defined  as 
positive  when  they  introduce  positive  force  contributions  along  the  X,  Y,  and  Z  axes,  respectively, 
at  small  angles  of  attack  and  sideslip  conditions.  Sign  convention  for  control  surface  deflections  are 
shown  in  Figures  6a-c. 

4.  AEROELASTIC  ANALYSIS 


The  aerodynamic  analysis,  like  the  structural  analysis,  is  based  on  a  finite  element  approach. 
The  finite  aerodynamic  elements  are  strips  or  boxes  on  which  there  are  aerodynamic  forces.  There 
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are  two  major  points  to  be  considered.  The  aerodynamic  elements,  even  for  rather  complex 
vehicles,  tend  to  be  regular  arrays.  Thus,  MSC/NASTRAN  aeroelastic  solution  generates  the 
arrays  of  aerodynamic  elements.  In  particular,  the  aerodynamic  elements  for  the  lattice  methods 
are  arrays  of  trapezoidal  boxes  whose  sides  are  parallel  to  the  airflow.  These  can  be  described 
simply  by  defining  properties  of  the  array  (panel).  Because  the  grid  points  defining  the  structural 
elements  usually  will  not  coincide  with  the  grid  points  defining  the  aerodynamic  elements,  provision 
has  been  made  to  generate  equations  for  interpolating  from  the  former  to  the  latter.  This  interpola¬ 
tion  is  a  key  feature  since  it  allows  the  choice  of  structural  and  aerodynamic  elements  to  be  based 
on  structural  and  aerodynamic  considerations  independently. 

Unsteady  aerodynamic  forces  are  generated  when  the  flow  is  disturbed  by  the  moving  (elastic) 
structure,  or  when  the  flow  itself  is  unsteady,  as  in  the  case  of  atmospheric  turbulence.  In  the  former 
case,  theory  leads  to  a  matrix  that  relates  the  forces  acting  on  the  structure  to  the  deflections  of  the 
structure.  State-of-the-art  methods  which  involve  interactions  among  aerodynamic  elements  are 
available  only  for  steady-state  sinusoidal  motion. 

4.1  Interconnection  of  the  Structure  with  Aerodynamics 

There  are  two  sets  of  displacements  for  the  analysis  of  static  aeroelasticity  in  MSC/NASTRAN. 
One  set  is  the  dependent  displacements  set,  Uk,  which  is  determined  at  a  set  of  points  whose 
location  is  determined  by  the  aerodynamic  theory.  The  other  set  is  the  independent  structural 
displacements  set,  u0,  which  consists  of  structural  grid  point  displacements  in  the  global  coordinate 
system.  The  dependent  displacements ,  uk,  are  interpolated  from  the  structural  displacement  ug  by 
the  following  relation. 

(Uk)  =  [Gkg]{Ug}  (1) 

In  order  to  obtain  the  interpolation  matrix  [G^],  three  types  of  splines  are  available  in 
MSC/NASTRAN.  These  are  out-of-plane  surfaces  spline,  out-of-plane  linear  spline,  and  in-plane 
linear  spline.  Out-of-plane  motion  consists  of  displacements  normal  to  the  plane,  and  rotations 
about  axes  parallel  to  the  plane.  In  this  study,  out-of-plane  displacements  (normal  displacement) 
over  the  wing  surface,  stabilator  surface,  and  vertical  fin  surface  are  interpolated  by  using  the  above 
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equation  (1).  Therefore,  the  out-of-plane  surface  spiine  is  used  for  this  analysis.  The  detailed 
explanation  of  a  surface  spline  method  is  given  in  [8]. 

4.2  Response  Solution 
The  structural  equation  is 


[K*]{ufl}  =  {F}  (2) 

where  [Kb]  is  the  structural  stiffness  matrix,  {Ug}  is  the  displacement  vector,  and  {F}  is  the  applied 
load  vector.  The  force  and  moments  are  represented  by 

{F}  =  {Fk}  +  {Ffl}  (3) 

where  {Fr}  is  the  applied  load  vector  that  contributes  to  rigid  body  motions  and  control  surfaces, 
and  {Fg}  is  the  additional  load  vector  due  to  the  change  from  the  original  deformation  (structural 
deformation).  The  aerodynamic  force  vectors  applied  at  rigid  body  motions  and  control  surfaces 
and  at  structural  grid  points  are  given  as 

{FkJ-qlSRjUAiir’lDjkKuk}  (4) 

{Fg}  =  qIGKg]T[Ski][Aiir1[Djk][Gkg]{ua}  (5) 

where  q  is  dynamic  pressure,  [Srq  is  a  force  transformation  matrix,  and  [Dkj]  is  a  displacement 
transformation  matrix.  The  force  vector  {Fr}  is  composed  of  two  loads, 

{Fr}  =  {Fs»}  +  {F8>  (6) 

where  {FM}  is  a  structural  load  vector,  and  {Fs}  is  a  static  aerodynamic  force  vector;  primarily,  the 
force  determined  by  the  planform  and  airfoil  configuration  of  the  wing  and  the  angle  of  attack. 

The  flexible  wing  geometry  can  be  expressed  in  terms  of  translational  and  rotational  deflection 
vectors, 


{U}  -  {Uk}  +  {Ug} 


(7) 
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where  (Uk)  is  the  undeformed  vector  at  the  aerodynamic  gird  points  and  (ug)  is  the  elastic 
deformation  vector  at  the  structural  grid  points. 


From  equation  (2)  through  (5),  we  obtain 

m  =  {Fk>  +  qIGkg]T [Skj][Ag]'1  [Dkj][Gkg]{Ug)  (8) 

Rearranging  equation  (8),  we  get 

[K]{Ug}  =  {Fk}  (9) 

where  the  system  stiffness  matrix  [K]  is  defined  as 

[K]»[K,]-[KaI  (10) 

and  the  aerodynamic  stiffness  matrix,  [Ka]  is  defined  as 

[Kal  =  qIGkg]T[Skj][Ajjr 1  [Dkj][Gkg]  (11) 

The  aeroelastic  response  solutions  are  obtained  from  equation  (9). 

5.  AERODYNAMIC  STABILITY  DERIVATIVES 

Stability  derivatives  are  aerodynamic  coefficients  nondimensionalized  by  reference  geometry 
appropriate  to  the  specified  derivative.  The  size  of  the  nondimensional  derivatives  is  normalized 
with  respect  to  aircraft  and  flight  conditions.  The  MSC/NASTRAN  aeroelastic  solution  provides  a 
complete  set  of  aerodynamic  estimates,  usually  as  a  function  of  angle  of  attack  for  the  various  flight 
conditions.  For  various  angles  of  attack,  the  downwash  equation  {W°j),  W2GJ  is  listed  on  the 
MSC/NASTRAN  DMI  card  option  (See  [9]  for  more  details).  In  the  longitudinal  case,  the  angle  of 
attack,  a,  and  stabilator  deflection,  5s,  are  the  two  incidence  variables, 


(12) 
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(13) 


where  S  and  c  are  the  reference  area  and  chord,  respectively.  The  factor  [rrv][M£  [introduces  the 
inertial  relief  effects  into  the  unrestrained  derivatives.  For  the  case  of  longitudinal  pitching  with  pitch 
rateq. 


(14) 


and 


KKFKHsfc}  (15) 

Similar  relationships  are  available  for  the  lateral-directional  stability  derivatives.  The  rigid 
stability  derivatives  are  appropriately  found  by  the  user’s  selection  of  a  small  value  of  q  and  do  not 
need  any  special  consideration. 


5.1  Longitudinal  Stability  Derivatives 


A  characteristic  of  fixed-wing  aircraft  longitudinal  stability  is  the  tendency  to  maintain  the  trim 
angle  of  attack.  Longitudinal  stability  derivatives  are  determined  by  the  aircraft  mass,  inertia 
aerodynamic,  and  geometric  layout  parameters.  The  most  important  derivatives  for  conventional 
aircraft  are:  Qa ,  Cm,.  Cm* ,  and  .  A  summary  of  these  longitudinal  terms  is  given  in  Table  2. 


Table  2.  Summary  of  Longitudinal  Terms 


Terms 

DescriDtion 

Affected  Bv 

Longitudinal  Static 
Stability 

Cma 

Tendency  of  aircraft  to  return 
to  trim  angle  of  attack 

Size,  location  of  stabilator, 
wing  pitching  moment  and 
fuselage  component  and 
relative  positions  of  c.g.  and 
center  of  lift  once 
configuration  is  fixed 

Pitch  Damping 
Coefficient 

Cmq 

A  measure  of  the  moment 
resistina  oitchina  motion 

Size  and  location  of  stabiiator 
and  air  density 

Angle  of  Attack 
Damping  Coefficient 

Qna 

A  measure  of  the  moment 
created  due  to  delay  of 
downwash  effect  on  tail  after 
wing  angle  of  attack  has  been 
chanced 

Size  and  location  of  stabiiator 
in  relation  to  wing 

Lift  Curve  Slope 

Lift  force  due  to  change  in 
angle  of  attack  (wing  or  tail) 

Area,  airfoil  shape,  aspect 
ratio  and  sweepback  angle  of 
wing  or  tail 

5.2  Lateral-Directional  Stability  Derivatives 

Lateral-directional  stability  is  generally  coupled  for  fixed-wing  aircraft  For  example,  while 
longitudinal  motion  is  restricted  to  a  single  axis  (pitch),  the  lateral-directional  motion  involves  two 
axes  (roll  and  yaw).  Also,  longitudinal  stability  derivatives  all  pertain  to  the  pitch  axis,  whereas 
certain  lateral  and  directional  derivatives  are  coupled  between  the  roll  and  yaw  axis.  Single 
derivative  variation  can  drastically  impact  the  lateral-directional  stability.  The  most  important 
lateral-directional  stability  derivatives  for  conventional  aircraft  are  C\p,  C^,  C^,  C^, ,  , 

and  Cyp.  A  summary  of  these  lateral-directional  terms  is  given  in  Table  3. 


Table  3.  Summary  of  Lateral-Directional  Terms 


Term 

OescriDtion 

Affected  Bv 

Directional  Static 
Stability  Coefficient 

Cnp 

Tendency  of  aircraft  to  align 
(weather  clock)  into  relative 
wind 

Primarily  vertical  tail  area  and 
distance  aft  of  c.g.  and  side 
area  distribution  force 

Yaw  Damping 
Coefficient 

Cry 

Yawing  moment  coefficient 
due  to  yaw  rate.  Created  by  (3 
changes  due  to  the  rotational 
velocity.  Acts  in  a  direction  to 
oppose  yaw  rate  (similar  in 
concept  to  pitch  damping 
coefficient  and  roll  damping 
coefficient) 

Size  and  placement  of  vertical 
tail  and  air  density 
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Yaw  due  to  Lateral 
Control  Coefficient 

c" * 

Yawing  moment  coefficient 
due  to  unbalanced  drag  of 
wings  with  aileron  deflection. 

Adverse  yaw  away  from 
desired  turn.  Pro  verse  yaw 
into  desired  turn 

Type  of  lateral  control  (aileron, 
etc)  and  control  size  and 
distance  from  c.g..  May 
affected  by  angle  of  attack  as 
airflow  over  vertical  tail  is 
blanked.  Effect  can  be 
minimized  with  aileron/rudder 
interconnect 

Yaw  due  to  Roll 
Rate  Coefficient 

CnP 

Yawing  moment  coefficient 
due  to  rotation  of  lift  vectors  of 
each  wind  durinq  roll 

Wingspan  and  lift  curve  slope 
(Cl  vs.  a) 

Rudder 

Effectiveness 

Coefficient 

Cnsr 

Yawing  moment  coefficient 
due  to  rudder  deflection 

Rudder  surface  area  and 
distance  of  rudder  from  c.g. 

Dihedral  Effect 

a 

Rolling  moment  coefficient 
due  to  sideslip.  Positive 
dihedral  -  roll  away  from 
sideslip.  Negative  dihedral  - 
roll  toward  sideslip  anale 

Geometric  dihedral  of  wings, 
location  of  wings  on  fuselage, 
and  effect  often  vary  with 
angle  of  attack 

’Roll  Damping 
Coefficient 

Rolling  moment  coefficient 
due  to  roll  rate.  Created  by  a 
changes  due  to  the  rotational 
velocity  "p“.  Acts  in  a  direction 
to  oppose  roll  rate  (similar  in 
concept  to  pitch  damping 
coefficient  and  yaw  damping 
coefficient) 

Lift  curve  slope  (CLvs.  a), 
wingspan  and  area,  and  air 
density 

Roll  due  to  Yaw 
Rate  Coefficient 

Cv 

Rolling  moment  coefficient 
created  by  one  wing  moving 
through  air  faster  than  other 
when  aircraft  vawed 

Lift  curve  slope  (Cl  vs.  a), 
and  wingspan 

2Roll  due  to  Rudder 
Deflection  Coefficient 

K9 

Rolling  moment  coefficient 
due  to  force  on  rudder  surface 
not  actino  throuqh  roll  axis 

Rudder  surface  area  and 
vertical  location  of  rudder 

Aileron 

Effectiveness 

Coefficient 

Cl* 

Rolling  moment  coefficient 
due  to  aileron  deflection 

Aileron  surface  area  and 
distance  from  ailerons  to  roll 
axis 

Side  Force 
Coefficient 

C„ 

Side-force  derivative 
coefficient  due  to  sideslip 

Fuselage  shape  and  aircraft 
geometric  layout 

NOTE:  ’Determines  roll  mode  time  constant  (how  quickly  roll  rate  reaches  steady  state) 
2Usually  a  minor  coupling  factor  on  most  aircraft 


5.3  Definition  of  Force  and  Moment  Coefficients 

The  principle  force  and  moment  coefficients  are  the:  lift  coefficient.  Cl;  rolling  moment 
coefficient  Q;  pitching  moment  coefficient  Cm;  yawing  moment  coefficient,  Cn;  and  the  normal  and 
axial  force  coefficients,  Cz  and  Cx.  The  longitudinal  derivatives  and  the  lateral  derivatives  are 
provided  as  the  output  as  shown  in  Tables  4-7.  Also,  the  MSC/NASTRAN  aeroelastic  solution  has 
the  capability  of  computing  the  aerodynamic  stability  derivatives  for  control  surface  deflection. 

Aerodynamic  coefficients  and  derivatives  are  used  to  calculate  the  aerodynamic  forces  and 
moments  acting  on  the  aircraft  Due  to  the  strong  dependence  of  the  and  derivatives  on 
the  angle  of  attack,  *hey  are  represented  as  angle  of  attack  dependent  coefficients  so  that: 

Cmi  *  Cmaflr 

CLi*CIaa 

and  the  lift  coefficient  is 

Cl  *  Cl,  +Cz^fi|  +•  C^q 

5.3.1  Baseline  Aircraft  Configuration 

Aerodynamic  force  and  moment  expansions  of  the  baseline  aircraft  configuration  result  in  the 
classical  description  of  the  symmetrical  aircraft  where  the  longitudinal  and  lateral  directional  forces 
and  moments  are  developed  independently.  The  baseline  aircraft  is  described  by  the  following  force 
and  moment  coefficients. 

Cy  “Cy^  +  Cy^^a  +Cy$.£  +CypP  +  Cyrr  (193) 

Cz  -  Cj,,  +  Cza<* + CZi8  fij  +  C^q 
C,  _  C,^+  Ct^  5a  +  $  +  QpP  +  C^r 

Cn  *  Cn^/J  +  ^  +  Cn$,  5-  +  CnpP  +  Cry  T 


(19b) 

(19c) 

(19d) 


(16) 

(17) 

(18) 
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5.3.2  Asymmetric  Wing  Damage  Configuration 

Aerodynamic  force  and  moment  expansions  for  the  asymmetric  wing  damage  configuration 
represent  a  significant  departure  from  the  baseline  aircraft.  This  configuration  is  characterized  by 
considerable  pitch/yaw  coupling  which  is  reflected  in  the  aerodynamic  model  which  is  presented 
below. 


Cy  »  Cyjj/3  +  Cy^  <^|  +  Cy^,  $  +  CypP  +  Cy,!  (203) 

C2  -  +  C  Zaa + CZpfi + Cz<5>  5,  +  Cz<5g  Sg  +  C^q  (20b) 

C,  «•  Cj  aa + C^fi + Q*  $1  +  C^,  £  +  CipP  +  C^q  +  C^r  (20c) 

Cn  *  Onpfi  +  Cn^  5a  +  Cn^,  Sr  +  C^p  +  C|yr  (20d) 


As  shown  in  the  above  equations,  to  calculate  the  rolling  moment  coefficient  of  asymmetric 
wing  damage  aircraft,  two  more  terms,  Ciaa  and  C^q ,  are  added  to  Eq.  (1 9c)  and  are  shown  in 
Eq.  (20c).  Two  terms,  Qa  and  ,  are  near  to  zero  for  undamage  wing  (baseline)  as  shown  in 
Tables  4-7.  So  we  can  neglect  these  two  terms  for  Eq.  (19c). 


6.  DISCUSSION  OF  RESULTS 

Results  are  obtained  for  five  aircraft  models,  one  baseline  (undamaged)  and  four  with 
asymmetric  wing  damages.  Four  subsonic  Mach  numbers  at  a  sea  level  flight  condition  are 
considered,  resultinc  in  the  dynamic  pressures  values  listed  in  Table  8. 

The  level  flight  condition  was  modeled  after  that  of  [9].  The  trim  solution  is  obtained  only  for 
subsonic  flight  in  quasi-steady  equilibrium.  In  this  study,  the  pitch  rate,  rudder  deflections,  pitch 
accelerations,  roll  accelerations,  yaw  accelerations,  lateral  accelerations,  and  vertical  accelera¬ 
tions,  are  all  set  to  zero.  Also,  to  compare  the  damage  case  4  (aileron  removed)  with  the  baseline 
case  and  the  other  three  damage  cases,  aileron  deflections  are  also  set  to  zero.  With  these  input 
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trim  conditions,  the  angle  of  attack,  sideslip  angle,  stabilator  rotation,  roll  rate,  and  yaw  rate  are 
determined. 

Table  9  shows  static  aeroelastic  effects  between  the  five  cases.  Results  are  obtained  for  a  1 0° 
angle  of  attack  and  a  Mach  number  of  0.7.  The  maximum  deflection  is  4.5168  in.  for  the  baseline 
wing  and  2.3947  in.  for  wing  damage  case  3.  The  roiling  moment  coefficient  is  smallest  for  the 
baseline  and  increases  with  wing  damage.  This  change  is  due  to  the  different  aerodynamic  load 
distributions  on  the  wing  surfaces.  The  lift  coefficient  calculation  is  based  on  equation  (18).  Lift 
coefficient  varies  linearly  with  wing  damage  due  to  the  loss  of  lifting  surface. 

Table  1 0  shows  trim  values  for  the  five  cases  at  Mach  0.7.  Here  the  roll  rate(p)  and  yaw  rate(r) 
increase  with  wing  damage.  It  is  shown  that  wing  damage  results  in  a  higher  roll  and  yaw  rate  than 
the  baseline  wing  for  trimmed  flight. 

The  effects  of  wing  damage  on  the  aerodynamic  coefficients  and  derivatives  are  shown  in 
Tables  4-7  at  Mach  numbers  0.4,  0.5,  0.6,  and  0.7  at  an  angle  of  zero  lift.  Among  the  variables 
shown  in  this  table,  Qa,  Qq,  and  are  the  most  significantly  affected.  This  is  largely  due  to 
longitudinal/lateral  coupling  effects.  These  terms  are  near  zero  for  the  baseline  aircraft. 

Table  1 1  shows  the  effects  of  normal  longitudinal  coupling  as  compared  with  lateral-directional 
motions.  Two  typical  derivatives,  Cma  and  ,  are  listed  in  this  table.  The  increments  of  weight 
loss  due  to  damaged  wings  and  their  locations  produce  asymmetric  inertial  terms  and  changes  the 
normal  axes  system  (the  center  of  gravity).  With  the  changes  in  aircraft  weight  and  inertia  due  to 
removal  of  wing  structure,  changes  in  aerodynamic  coefficient  derivatives  are  expected. 

Figure  7  shows  the  aileron  effectiveness.  The  rolling  moment  due  to  aileron  deflection  is 
significantly  affected  by  aileron  area  removal  (or  distance  from  aileron  to  roll  axis).  For  damage 
cases  3  and  4,  the  entire  right  aileron  is  removed. 

Figure  8  shows  that  wing  damage  causes  a  reduction  in  lift  slope  with  very  little  change  in  the 
angle  of  zero  lift.  Although  the  wing  is  considered  as  elastic,  the  lift  coefficient  is  linear  for  this  angle 
of  attack  range. 
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Figure  9  shows  the  changes  in  the  lift  curve  slope  with  respect  to  the  angle  of  attack,  a,  for  the 
baseline  and  the  four  damaged  wing  configurations,  it  is  shown  that  this  term  increases  with  Mach 
number. 

Comparison  of  wing  damage  data  in  Figures  lOa-b,  11a-b,  and  12a-b  shows  large  systematic 
changes  in  the  rolling  moment  (Q)  which  are  dependent  mainly  on  angle  of  attack  a.  At  a=  5.0°, 
Figure  10b  shows  increased  rolling  moments  as  outboard  portions  of  the  wing  are  removed. 
Damage  case  3  (30.3  %  right  wing  off)  has  the  largest  rolling  moment  among  the  four  cases.  This 
trend  was  not  observed  for  small  a  as  shown  in  Figure  10a  (o=0.4°).  While  yav  i  moment  and 
side  force  are  not  affected  by  angle  of  attack  for  the  damage  wing  configurations. 

The  yawing  moment  coefficient  Cn,  shown  in  Figures  lla-b,  is  of  particular  interest:  the  side 
force  coefficient  Cy,  shown  in  Figures  12a-b,  changes  very  little,  and  the  rolling  moment  coefficients 
behave  as  would  have  been  expected;  that  is,  increasing  lift  with  the  lateral  center  of  pressure 
moved  to  the  left  of  the  origin  would  be  expected  to  increase  the  rolling  moment.  The  quite  regular 
reduction  in  the  lift  coefficient  slope  with  respect  to  the  angle  of  attack  for  the  baseline  aircraft  and 
the  four  wing  damage  configurations  are  summarized  in  Figure  9. 

Figures  13a-b  show  that  with  sideslip,  there  is  very  little  influence  on  pitching  moment 
coefficient  Figure  13b  also  shows  that  (at  Mach  0.7  and  a=5.0°)  wing  damage  significantly 
influenced  the  pitching  moment. 

Wing  damage  reduces  both  the  wing  area  and  the  lateral  distance  from  the  aircraft  centerline 
to  the  wing  center  of  pressure.  The  rolling  moment  is  calculated  as  the  difference  in  the  products 
of  the  lift  and  distance  to  the  spanwise  center  of  pressure  for  the  damaged  and  undamaged  wings. 
The  calculated  rolling  moment  coefficient  for  the  baseline  and  the  four  damaged-wing  aircraft  is 
plotted  in  Figure  14  at  Mach  =  0.4, 0.5, 0.6,  and  0.7,  and  at  an  angle  of  attack  of  10.0°. 

The  most  important  effect  of  wing  damage  is  shown  in  Figures  15-17,  where  the  rolling  moment 
coefficient  is  seen  to  be  strongly  affected  by  angle  of  attack.  The  rolling  moment  coefficients  (Q), 
pitching  moment  coefficients  (Cm),  and  yawing  moment  coefficients  (Cn)  are  plotted  vs.  sideslip 
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angle  for  three  different  angles  of  attack  at  M=0.7  for  the  damage  case  3  (30.3  %  right  wing  off)  and 
the  baseline  aircraft.  At  angles  of  attack  of  0.4°,  5.0°,  and  10.0°,  the  baseline  aircraft  rolling  moment 
is  not  changed.  The  three  baseline  curves  are  all  at  the  same  location.  However,  the  rolling  moment 
with  damage  increases  with  angle  of  attack.  It  is  shown  that  the  slope  of  the  rolling  moment 
coefficient  vs.  sideslip  (c^)  is  not  greatly  affected  by  wing  damage  or  angle  of  attack  even  though 
the  magnitude  of  the  induced  rolling  moment  increases  significantly.  For  the  baseline  aircraft  the 
C|  vs.  0  variation  is  shown  to  be  insensitive  to  angle  of  attack  from  zero  to  10°  with  a  value  of 
approximately  -0.096112. 

Figure  16  shows  that  the  pitching  moment  (Cm)  is  reduced  with  angle  of  attack  and  wing 
damage  due  to  reduction  in  total  lift  and  shift  of  the  center  pressure.  Also,  the  pitching  moment  is 
influenced  very  little  by  the  amount  of  sideslip.  Figure  1 7  shows  that  the  yawing  moment  coefficient 
(Cn)  is  strongly  affected  by  sideslip,  but  not  by  wing  damage. 

Figures  1 8a  and  1 8b  show  the  damping-in-roll  as  an  elastic  wing  and  rigid  wing  for  the  baseline 
and  damage  cases  at  Mach  numbers  of  0.4, 0.5, 0.6,  and  0.7.  Flexibility  effects  increase  with  Mach 
number  for  all  cases. 
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Table  4.  Aerodynamic  Coefficient  Derivatives  at  Mach  0.4 


Mach  =  0.4 

Model 

Baseline 

Damaqe  1 

Damaae  2 

Damaae  3 

Damaae  4 

4.21504E+00 

3.79706E+00 

3.52388E+00 

2.91184E+00 

4.04093E+00 

-8.94627E+00 

-7.57379E+00 

-9.06388E+00 

-7.84992E-06 

4.17877E-02 

-7.39723E-02 

-2.58857E-01 

-2.73841  E-01 

(291 

5.71944E-01 

5.69672E-01 

5.67432E-01 

5.62092E-01 

5.70935E-01 

\smm 

-1.44600E+00 

-1.22174E+00 

-1 .08762E+OU 

-1. 00846 E+00 

-1.35119E+00 

M 

-7. 1071  IE-04 

1.41237E-01 

1.98113E-01 

2.28794E-01 

3.58894E-02 

H 

-1.01680E+01 

HSHSIBI 

-9.5841 9E+00 

2.38867E-04 

2.72647E-02 

-4.23605E-02 

-1.39133E-01 

-1.444948E-01 

(231 

1.37514E+00 

1.37474E+00 

1.37356E+00 

1.37260E+00 

1.37545E+00 

Cla 

1. 67761 E-04 

1.26804E-01 

1.89476E-01 

2.82202E-01 

3.57967E-02 

-7.81310E-02 

-7.74624E-02 

-7.851 18E-02 

-7.851 76E-02 

-7.85873E-02 

m 

-3.61364E-01 

-2.83867E-01 

-2.52809E-01 

-2.21798E-01 

-3.48072E-01 

<*o 

7.401 25E-04 

1.94531  E-01 

2.901 00E-01 

4.36061  E-01 

1.04029E-01 

Cv 

9.74985E-02 

9.6631  IE-02 

9.71206E-02 

9.7371 4E-02 

9.80339E-02 

M 

1.10443E-01 

1.18387E-01 

8.84361  E-02 

5.67570E-02 

5.59648E-02 

H 

-1.56579E-04 

-1.12863E-03 

-1.59978E-03 

-2.34888E-03 

-5.31615E-04 

4.56876E-02 

4.51982E-02 

4.51368E-02 

4.52327E-02 

4.58502E-02 

BHM 

-4.11054E-01 

-4.10737E-01 

Btt 

-1 .08371 E-01 

-1.15468E-01 

-1.20547E-01 

-1.29159E-01 

-1.11145E-01 

E2NMK 

2.07602E-01 

2.07735E-01 

2.07793E-01 

2.0791  IE-01 

2.07636E-01 

2.3331 7E-05 

5.43758E-03 

9.58457E-03 

2.04680E-02 

3.14829E-03 

M 

2.0031  IE-01 

2.00420E-01 

2.00366E-01 

2.00483E-01 

2.00323E-01 

eeshp 

5.44556E-02 

5.79462E-02 

6.03200E-02 

6.42365E-02 

5.57621  E-02 

EH  1 

-2.68094E-04 

7.97308E-03 

1.42244E-02 

3.10524E-02 

9.07929E-03 

-2.5921 3E-01 

-2.59351  E-01 

-2.59358E-01 

-2.59479E-01 

-2.59231  E-01 

n 

9.9501 7E-03 

1.06599E-02 

8.17248E-03 

4.77960E-03 

5.01351  E-03 

H 

6.86758E-05 

1. 52971 E-05 

-2.241 33E-05 

-1.33222E-04 

2.33446E-05 

rail 

-1.22296E-01 

-1.22362E-01 

-1.22390E-01 

-1.22447E-01 

-1.2231  IE-01 

Tab!?  5.  Aerodynamic  Coefficient  Derivatives  at  Mach  0.5 


Mi 

ach  =  0.5 

Model 

Baseline 

Damaae  1 

Damaae  2 

Damaae  3 

Damage  4 

SB 

4.57659E+00 

4.10584E+00 

3.81836E+00 

3.17836E+00 

4.37825E+00 

n 

-1.16997E+01 

-1.09620E+01 

-1.05282E+01 

-9.4941 2E+00 

-1.11010E+01 

-1. 45461 E-04 

4.76725E-02 

-2.79884E-01 

-2.88440E-01 

n 

1. 06051 E+00 

1.04891  E+00 

1.04221  E+00 

1.02540+00 

1.05214E+00 

Cma 

-1.77552E+00 

-1.50672E+00 

-1.35713E+00 

-1.27100E+00 

-1.65597E+00 

Cn* 

-2.85686E-04 

1.66015E-01 

2.34028E-01 

2.63505E-01 

4.18294E-02 

-1.31054E+01 

-1.26693E+01 

-1.24240E+01 

-1.21474E+01 

1.27511E+01 

&39H 

4.8371 4E-04 

3.36043E-02 

-1.62964E-01 

-1.64722E-01 

(■S 

2.10135E+00 

2.09495E+00 

2.091 33E+00 

2.08780E+00 

2.09689E00 

mm 

-1.96970E-04 

1.35930E-01 

2.05927E-01 

3.001 17E-01 

3.91342E-02 

-7.73243E-02 

-7.67339E-02 

-7.59866E-02 

-7.64892E-02 

-7.79386E-02 

^ _ 

-3.81447E-01 

-2.96990E-01 

-2.64953 E-01 

-2.35036E-01 

-3.67553E-01 

EBB 

-5.36550E-04 

2.07862E-01 

3.19944E-01 

4.75023E-01 

1.14177E-01 

Cv _ 

9.69865E-02 

9.57009E-02 

9.52461  E-02 

9.5831 9E-02 

9.771 53E-02 

smm 

1.12176E-01 

1.19102E-01 

9.15306E-02 

5.90544E-02 

5.701 52E-02 

Cu. _ 

1.05286E-04 

-2.98373E-03 

-4.62703 E-03 

-1.35028E-03 

EBB 

4.75004E-02 

4.6621 3E-02 

4.64843E-02 

4.66492E-02 

4.77401  E-02 

BB 

-4.1 6073 E-01 

EHSEHEfll 

-4.15912E-01 

EBB 

-1.06301  E-01 

-1.541 47E-01 

-1.19302E-01 

-1.27519E-01 

-1.09260E-01 

mm 

2.11482E-01 

2.11673E-01 

2.11759E-01 

2.111908E-01 

2.11 521  E-01 

2.05693E-06 

5.61868E-03 

1.04268E-02 

2.07636E-02 

3.391 01  E-03 

b 

Iwl 

2.02581  E-01 

2.02670E-01 

2.02802E-01 

2.02496E-01 

EBB 

5.36605E-02 

5.75998E-02 

5.99385E-02 

6.36624E-02 

5.50501  E-02 

EBB 

bebbhsi 

7.86724E-03 

1.57873E-02 

3.30885E-02 

9.84587E-03 

C"t _ 

-2.6261 6E-01 

-2.62779E-01 

-2.62887E-01 

-2.63021  E-01 

-2.62633E-01 

M 

9.94729E-03 

1.03764E-02 

8.64933E-03 

5.04588E-03 

5.01822E-03 

EB 

1.02969E-04 

-1.80294E-04 

-4.92736E-04 

-5.13809E-05 

SB 

-1.25256 E-01 

-1.2535  IE-01 

-1.25394E-01 

-1.25469E-01 

-1.25273E-01 
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Table  6.  Aerodynamic  Coefficient  Derivatives  at  Mach  0.6 


Mach  =  0.6 


Model 

Baseline 

Damaqe  1 

Damaqe  2 

Damaqe  3 

Damaqe  4 

mm 

4.69352E+00 

4.19492E+00 

3.89357E+00 

3.32482E+00 

4.501 14E+00 

mm 

-1.13998E+01 

-1.06459E+01 

-1.02122E+01 

-9.321 84E+01 

-1.08420E+01 

mm 

3.58901  E-04 

4.79920E-02 

H5S5SI 

-2.59486E-01 

-2.64279E-01 

mm 

8.53428E-04 

8.44764E-01 

8.40099E-01 

8.30662E-01 

8.48295E-01 

Cnv* 

-1.61352E+00 

-1 .34461 E+00 

-1.33607E+00 

-1.19597E+00 

-1.53106E+00 

Cmp 

-1. 57721 E-03 

1.72766E-01 

2.08004E-01 

2.36287E-01 

2.78869E-02 

B5H 

-1 .21637E+01 

-1.17359E+01 

-1.15910E+01 

-1.14989E+01 

-1.19522E+01 

BH 

8.88345E-04 

3.36563E-02 

esssssi 

-1.08155E-01 

-1.05554E-01 

m 

1.72049E+00 

1.71600E+00 

1.71501  E+00 

1.71362E+00 

1.71907E+00 

1.46705E-01 

2.22338E-01 

3.12660E-01 

4.0921 0E-02 

SHI 

-7.91327E-02 

-7.9681 6E-02 

Clo 

-4.04776E-01 

-3.11563E-01 

-2.77903 E-01 

-2.50060E-01 

-3.90565E-01 

On 

1.81228E-03 

2.17391  E-01 

3.35373E-01 

4.81226E-01 

1.18673E-01 

Cfe _ 

1.02180E-01 

9.95691  E-02 

9.88831  E-01 

9.95305E-02 

1.03046E-01 

Mi 

1.19518E-01 

1.26233E-01 

9.97864E-02 

6.44079E-02 

6.37744E-02 

mm 

BUSSES! 

-3.17915E-03 

-4.31597E-03 

-5.58684E-03 

-1.79059E-03 

mm 

5.031 90E-02 

4.88097E-02 

4.85593E-02 

4.87702E-02 

5.06289E-02 

mm 

-4.16945E-01 

EHH3351 

KSESSSl 

-4.17038E-01 

Ml 

-1.06883E-01 

-1.15842E-01 

-1.20564E-01 

-1.27375E-01 

-1.09686E-01 

mu 

2.16641 E-01 

2.16942E-01 

2.17074E-01 

2.17316E-01 

2.16715E-01 

mm 

-2.661 25E-04 

4.95746E-03 

9.9421 3E-03 

1.62459E-02 

2.4271 7E-03 

EP  1 

2.02555E-01 

2.02757E-01 

2.02899E-01 

2.03094E-01 

2.02597E-01 

Crv, _ 

5.41 51  IE-02 

5.83790E-02 

6.0661 7E-02 

6.35370E-02 

5.54097E-02 

3.12841  E-04 

7.391 30E-03 

1.57035E-02 

2.60384E-02 

7.94506E-03 

mm 

-2.641 87E-01 

WSBSSSA 

-2.64621  E-01 

-2.64867E-01 

-2.64243E-01 

7.57369E-03 

7.64696E-03 

6.27860E-03 

3.67670E-03 

3.78637E-03 

mm 

-2.66566E-04 

-3.7361 7E-04 

-4.72809E-04 

B3S5SB1 

-3.4921 2E-04 

mm 

-1.29370E-01 

-1.29519E-01 

-1.29585 E-01 

-1.29709E-01 

-1.29403E-01 
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Table  7.  Aerodynamic  Coefficient  Derivatives  at  Mach  0.7 


Mach  =  0.7 


Model 

Baseline 

Damaqe  1 

Damaqe  2 

Damaqe  3 

Damaqe  4 

Cz« _ 

5.15479E+00 

4.58625E+00 

4.25004E+00 

3.50202E+00 

4.87709E+00 

HI 

-1.11120E+01 

-1.03106E+01 

-1.02645E+01 

O  1 

-2.2221 6E-04 

4.2901 7E-02 

TBSBB&M 

-3.46944E-01 

-3.49901  E-01 

H 

7.49472E-01 

7.39272E-01 

7.33980E-01 

7.18922E-01 

7.40505E-01 

Crria 

-1.98009E+00 

-1.62726E+00 

-1 .46588E+00 

-9.14994E-01 

-1.73563E+00 

Or* 

3.25451  E-04 

2.11031E-01 

3.03260E-01 

4.42684E-01 

8.89820E-02 

-7.92257E+00 

-8.97036E+00 

H 

4.31118E-02 

-2.98399E-01 

-3.03681  E-01 

E3HH 

1.31379E+00 

1.30749E+00 

1.30510E+00 

1.29420E+00 

1.30552E+00 

-1 

1.07583E-03 

1.76055E-01 

2.58507E-01 

4.06247E-01 

5.98702E-02 

mi 

-9.21453E-02 

-9.15835E-02 

-9.68908E-02 

Hi 

-4.49736E-01 

-3.42647E-01 

-3.03799E-01 

-2.65389E-01 

-4.30996E-01 

H 

-1.89509E-04 

2.46285E-01 

3.63092E-01 

6.04267E-01 

1.70703E-01 

Cy  1 

1.23854E-01 

1.18732E-01 

1.17994E-01 

1.19517E-01 

1.24789E-01 

m 

1. 39481 E-01 

1.49958E-01 

1.18647E-01 

7.38329E-02 

7.15809E-02 

SHU 

5.40273E-04 

-2.22091  E-03 

-3.36939E-03 

-5.52093E-03 

H 

6.15065E-02 

5.89076E-02 

5.8491 9E-02 

5.90831  E-02 

6.18118E-02 

H 

-4.17418E-01 

-4.1 8441  E-01 

-4.18878E-01 

-4.17640E-01 

H 

-1.02890E-01 

-1.1 3041  E-01 

-1 .1 91 66E-01 

-1.29941  E-01 

-1.07278E-01 

HI 

2.22029E-01 

2.22479E-01 

2.22683E-01 

2.22993E-01 

2.22203E-01 

H 

2.23320E-04 

7.73506E-03 

1.26650E-02 

3.07653E-02 

5.67604E-03 

BP  1 

2.01431  E-01 

2.01765E-01 

2.01919E-01 

2.021 16E-01 

2.01542E-01 

^Qa _ 1 

5.16275E-02 

5.65227E-02 

5.951 66E-02 

6.51267E-02 

5.38698E-02 

SHI 

-1.71128E-04 

1.06530E-02 

1.71238E-02 

4.71075E-02 

1.73886E-02 

SB 

-2.6381 6E-01 

-2.64247E-01 

-2.64445E-01 

-2.64706 E-01 

-2.63965E-01 

S!  1 

1. 39791 E-02 

1.47233E-02 

1.26119E-02 

7.06278E-03 

7.04896E-03 

mm 

1.05204E-04 

-4.6681 0E-04 

-1 .30435E-04 

EB55351 

-9.261 75E-05 

m 

-1.33899E-01 

-1.34113E-01 

-1.3421  IE-01 

-1.34354E-01 

-1.33987E-01 
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Table  8.  Dynamic  Pressure  Variables  -  Standard  Sea  Level 


Mach  Is 

lumber 

0.4 

0.5 

0.6 

0.7 

Dynamic  Pressure,  psi 

1.7811 

2.7821 

4.0075 

5.4571 

Table  9.  Comparison  of  Aeroelastic  Effects  for  Baseline  and  Damaged  Wings 


Mach  =  0.7 

Model 

BUHffS 

Damaae  2 

Damaae  4 

0.0  % 

8.0% 

13.73  % 

30.3  % 

6.63% 

24319.1 

23353.44 

22649.67 

20634.76 

23513.48 

Maximum  Deflection,  in 

4.5168 

3.5429 

3.1171 

2.3947 

4.2697 

Lift  coefficient 

0.8992 

0.8001 

0.7414 

0.6109 

0.8508 

Total  Lift,  lb 

119,334.4 

101.966.4 

91.638.2 

68.790.9 

109.170.8 

Anale  of  Attack,  dea 

10 

10 

10 

10 

10 

Rolling  Moment  Coefficient 

0.000187 

0.030712 

0.045096 

0.070868 

0.010444 

Table  10.  Trim  Values  at  Level  Flight  at  Mach  0.7 


Trim  at  Level  Fli 

ight  at  Mach  0.7 

BBSISSHB 

Damaae  2 

WSB&M 

_ 8a.  rad _ 

0.0 

0.0 

0.0 

HEHS9HI 

8.05331  E-05 

7.04628E-05 

6.3591 2E-05 

1.26446E-05 

7.31 528E-05 

8r.rad 

0.0 

0.0 

0.0 

0.0 

0.0 

a.  rad 

5.60787E-05 

6.351 18E-05 

7.24653E-05 

1.16160E-04 

5.79581  E-05 

_ B.  rad _ 

1.66542E-06 

7.37742E-07 

-2.70060E-06 

1.89746E-05 

3.44328E-06 

2.44259E-07 

3.3731 7E-05 

6.53781  E-05 

2.00482E-04 

8.04404E-06 

0.0 

0.0 

mmmm 

0.0 

0.0 

r,  rad/sec 

1.40055E-06 

9.53687E-05 

1.59679E-05 

7.71279E-05 

5.46505E-06 
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Figure  8.  Lift  coefficient  at  Mach=0.7.  a*&0^Q0M°.  6.0°.  8.0°.  and  1Q.Qt 
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Figure  9.  Lift  curve  slope  vs  Mach=Q.4.  Q.5.  Q.6.  and  Q.7. 
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SIDESLIP  ANGLE  PEG) 

Figure  11a.  Yawing  moment  coefficient  at  Mach=0.7. 


Figure  11b. 
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Figure  12a.  Side  force  coefficient  at  Mach=0.7.  a=C 
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Figure  12b. 
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» 13a.  Pitching  moment  coefficient  at  Mach=0.7 
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Figure  18a.  Damping-in-roll  elastic  wing. 


Figure  18b.  Damping-in-roll  rigid  wing. 
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7.  CONCLUSIONS 


This  study  analyzed  the  static  aeroelastic  response  of  an  aircraft  with  asymmetric  wing 
planforms  representative  of  combat  damage.  Trim  variables,  static  deformations,  air  loads,  and  the 
aerodynamic  coefficients  and  derivatives  are  obtained  for  both  undamaged  and  damaged  wing 
cases. 

The  following  conclusions  are  drawn  from  this  study. 

1 )  The  effect  of  longitudinal/lateral  coupling  for  an  asymmetric  wing  configuration  is  quite  significant 
with  wing  damage. 

2)  The  rolling  moment  coefficient  is  the  most  important  stability  derivative  in  the  study  of  aircraft 
with  wing  damage.  As  discussed  above,  the  effects  are  significantly  influenced  by  the  angle  of 
attack. 

3)  Based  on  trim  considerations,  possible  flight  attitudes  (angle  of  attack  and  sideslip)  were  found 
that  will  permit  straight  and  level  flight  for  a  contemporary  aircraft  configuration  with  a  major  loss  of 
the  aerodynamic  lifting  and  control  surfaces.  However,  unstable  dynamic  modes  may  exist  which 
could  render  the  aircraft  difficult  or  impossible  to  control.  To  address  this  possibility,  the  present 
results  will  be  used  with  an  aircraft  simulation  code  (ACSIM)  based  on  a  coupled  six  degree-of- 
freedom  aircraft  equation  of  motion  analysis. 

Wing-body-tail  interference  effects  are  being  considered  for  future  study  to  predict  the 
aerodynamic  effects  of  combined  damage  configurations,  i.e.,  fuselage  damage  combined  with  a 
damaged  wing,  a  damaged  stabilator,  and/or  a  damaged  vertical  tail. 


8.  RECOMMENDATIONS 


Through  this  research,  three  recommendations  for  future  study  are  suggested,  a)  Investigate 
the  effects  of  the  aeroelastic  response  of  a  damaged  stabilator  or  vertical  wing,  and  develop  a 
structural  and  aerodynamic  model  of  the  tail  group,  b)  Based  on  the  author’s  experience  with  finite 
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element  wing  models  (See  [5]  and  [6]),  the  box-beam  type  wing  model  should  be  optimized  as  a 
plate-beam  type  wing  model  to  save  computing  time  and  to  make  it  easier  to  build  a  finite  element 
model  for  damaged  wing  cases  for  static  aeroeiastic,  flutter,  and  dynamic  response  analysis  using 
the  MSC/NASTRAN  aeroeiastic  solution,  c)  A  wind  tunnel  test  should  be  planned  and  conducted 
to  develop  criteria  to  correlate  the  calculation  method  of  the  damaged  wing  aerodynamics. 
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APPENDIX  A: 

Data  Input  and  Generation  for  MSC/NASTRAN  Aeroelastidty 
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A1 .  Data  Input  and  Generation 


Aerodynamic  calculations  are  performed  using  a  Cartesian  coordinate  system.  By  user 
convention,  the  flow  is  in  the  positive  x-direction  and  x-axis  of  every  aerodynamic  element  and  must 
be  parallel  to  the  flow  in  its  undeformed  position.  The  basic  structural  coordinate  system  may  be 
defined  independently,  since  the  use  of  the  same  system  for  both  would  place  an  undesirable 
restriction  upon  the  description  of  the  structural  model;  and  any  MSC/NASTRAN  Cartesian  system 
may  be  specified,  as  long  as  the  flow  is  defined  in  the  direction  of  the  x-axis.  All  aerodynamic 
element  and  grid  point  data,  initially  defined  in  the  basic  coordinate  system,  will  be  transformed  to 
the  aerodynamic  coordinate  system.  All  the  global  (displacement)  coordinate  systems  of  the 
aerodynamic  grid  points  will  have  their  x-directions  in  the  flow  direction.  Their  z-directions  will  be 
normal  to  the  element  in  the  case  of  boxes,  and  parallel  to  the  aerodynamic  z-direction  in  the  case 
of  bodies. 

The  aerodynamic  grid  points  are  physically  located  at  the  center  of  the  boxes  for  the  lifting 
surface  theories,  and  at  the  centers  of  body  elements  for  the  Doublet-Lattice  method.  Permanent 
constraints  are  generated  for  the  unused  degrees  of  freedom.  A  second  set  of  grid  points  is  used 
only  for  the  element  identification  number.  For  any  panel,  the  box  numbers  start  with  the  panel 
identification  number  and  increase  consecutively. 

A2.  Aerodynamic  Modeling  in  MSC/NASTRAN 


Based  on  the  guidelines  explained  above  in  section  A1 ,  the  descriptions  of  input  data  to  create 
the  aerodynamic  model  are  shown  in  card  format  in  Tables  A1-A3.  These  all  option  cards  are  used 
in  the  BULK  DATA  CARD  [Ref.  9]. 


Table  A1 .  Aerodynamic  Elements 


Description 

AEFACT 

Specifies  lists  of  real  numbers  for  aeroelastic  analysis 

m 

Defines  relationship  between  AESTAT  and  AESURF  entries 
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AELIST 


AESTAT 


AESURF 


CAER01 


CAER02 


PAER01 


PAER02 


AEROS 


Defines  a  list  of  aerodynamic  elements  to  undergo  the  motion  prescribed  with 
the  AESURF  Bulk  Data  entry  for  static  aeroelasticity 


Specifies  rigid  body  motions  to  be  used  as  trim  variables  in  static  aeroelastici 


Specifies  an  aerodynamic  control  surface 


Defines  an  aerodynamic  macro  element  (panel)  in  terms  of  two  leading  edge 
locations  and  side  chords.  This  is  used  for  Doublet-Lattice  theory  for  subsonic 

aerodynamics 


Defines  aerodynamic  slender  body  and  interface  elements  for  Doublet-Lattice 

aerodynamics 


Defines  aerodynamic  panel  properties  in  the  Doublet-Lattice  method 


Defines  the  cross-sectional  properties  of  aer 


Basic  physical  data  for  static  aeroelastici 


Table  A2.  Aeroelastic  Response  Analysis 


Description 

TRIM 

Specifies  constraints  for  aeroelastic  trim  variables 

Table  A3.  Aerodynamic  to  Structure  Interconnection 


r 

Description 

SET1 

Defines  a  set  of  structural  arid  points  by  a  list 

SPLINE1 

Defines  a  surface  spline  for  interpolating  out-of-plane  motion  for  aeroelastic 

problems 

SPLINE2 

Defines  a  beam  spline  for  interpolating  panels  and  bodies  for  aeroelastic  - 

problems 
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1NV3  WWl- 


Aileron  Hinge  Line 


NOTE:  The  stabilator  is  an  all  moving  tail 


Rudder  Hinge  Line 


Rudder  Aero  Box 


Cl.  Introduction 


As  described  earlier  in  this  report,  the  static  aeroelastic  analysis  is  designed  to  obtain  both 
structural  and  aerodynamic  data.  The  structural  data  of  interest  include  loads,  deflections,  and 
stresses.  The  aerodynamic  data  include  stability  and  control  derivatives  and  trim  values.  The 
analysis  presupposes  a  structural  model,  an  aerodynamic  model,  and  the  interpolation  between 
the  two.  In  this  report,  MSC/NASTRAN  V67  on  the  ail  Cray  computer  is  used  to  analyze  this  study 
using  Sol  21  or  144  (for  super-element)  on  the  Executive  Control  Card  Deck. 

C2.  Option  for  the  Case  Control  Card 

The  static  deflections,  stress,  strain,  loads,  aerodynamic  pressures,  and  aerodynamic  forces 
are  obtained  as  part  of  the  options  shown  in  Table  Cl .  These  derived  quantities  of  interest  must  all 
be  requested  in  the  Case  Control  Deck.  The  various  flight  conditions  are  specified  for  each 
SUBCASE  on  TRIM  cards. 


Table  Cl .  Option  Card  for  Output 


Card 

DescriDtion  of  Output 

APRES 

Request  the  aerodynamic  pressures 

AEROF 

Reauest  the  aerodynamic  forces 

FORCE 

Request  the  structural  elements  loads 

DISP 

Reauest  the  structural  deflections 

STRESS 

Request  the  structural  element  stresses 

C3.  Option  for  the  Bulk  Data  Card 

In  the  Bulk  Data  Card  Deck,  PARAM  POST  is  used  to  create  an  external  database  that  can  be 
used  for  plotting.  The  POST  =  0  option  is  for  MSC/XL.  The  database  can  be  used  to  plot  analysis 
results;  for  example,  to  generate  contour  plots  of  stress,  strain,  and  displacements  on  an  aircraft 
wing  due  to  airloading. 
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The  stability  derivatives  are  obtained  as  part  of  the  solution  process  and  are  always  printed 
follows  request.  These  are  shown  in  Table  C2. 


Table  C2.  Option  Card  for  Stability  Derivatives 


Label 

DescriDtion 

ANGLEA 

Anqle  of  attack 

SIDES 

Anqle  of  sideslip 

PITCH 

Pitch  rate 

ROLL 

Roll  rate 

YAW 

Yaw  rate 

AILERON 

Deflection  of  aileron 

ELEV 

Deflection  of  elevator 

RUDDER 

Deflection  of  rudder 

URDD2 

Lateral  acceleration 

URDD3 

Vertical  acceleration 

URDD4 

Roll  acceleration 

URDD5 

Pitch  acceleration 

URDD6 

Yaw  acceleration 
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